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The devil’s staircase is a fractal structure that characterizes the ground state of one-dimensional 
classical lattice gases with long-range repulsive convex interactions. Its plateaus mark regions of 
stability for specihc hlling fractions which are controlled by a chemical potential. Typically such 
staircase has an explicit particle-hole symmetry, i.e., the staircase at more than half-hlling can 
be trivially extracted from the one at less than half hlling by exchanging the roles of holes and 
particles. Here we introduce a quantum spin chain with competing short-range attractive and 
long-range repulsive interactions, i.e. a non-convex potential. In the classical limit the ground state 
features generalized Wigner crystals that — depending on the hlling fraction — are either composed 
of dimer particles or dimer holes which results in an emergent complete devil’s staircase without 
explicit particle-hole symmetry of the underlying microscopic model. In our system the particle-hole 
symmetry is lifted due to the fact that the staircase is controlled through a two-body interaction 
rather than a one-body chemical potential. The introduction of quantum huctuations through a 
transverse held melts the staircase and ultimately makes the system enter a paramagnetic phase. 
For intermediate transverse held strengths, however, we identify a region, where the density-density 
correlations suggest the emergence of quasi long-range order. We discuss how this physics can be 
explored with Rydberg-dressed atoms held in a lattice. 


Introduction .—Systems with long-range interactions 
can feature intriguing physics that is not necessarily 
present in their short-range counterparts. For example, 
classical particles in a one-dimensional (ID) lattice in¬ 
teracting via repulsive inhnite-range convex potentials 
lower their interaction energy by assuming a distribution 
in space as uniform as possible. Remarkably this prop¬ 
erty makes the ground state conhguration independent 
of the actual details of the interactions, e.g. the specihc 
power-law of the interaction potential, and leads to the 
formation of a so-called generalized Wigner crystal [1, 2]. 
Furthermore, the permitted hlling fractions of the ground 
state conhguration versus the chemical potential form a 
fractal curve known as the complete devil’s staircase [3]. 

Recent years have seen a great success in emulating 
many-body systems with long-range interactions using 
ultracold gases. Here, crystalline structures which were 
originally studied in the context of solid state physics 
have shown to be present also in ensembles of cold 
trapped ions [4], polar molecules [5-7] and gases of Ryd¬ 
berg atoms [8-11]. In particular Rydberg gases have wit¬ 
nessed a recent experimental breakthrough in which the 
adiabatic preparation of a crystalline state and the onset 
of a staircase structure were shown [12]. A case currently 
much less studied is when long-range interactions feature 
competing attraction and repulsion. Systems with such 
interactions are known to exhibit intricate behavior and 
in fact appear in a number of contexts, e.g., Langmuir 
monolayers [13], classical huids [14, 15], mixtures of col¬ 
loids and polymers [16, 17] and ferrohuids [18]. 

In this work we explore such scenario in a cold atomic 
lattice gas in which interatomic interactions are con¬ 
trolled via a so-called double Rydberg dressing scheme 
[19-21]. This permits the realization of competing 


nearest-neighbor (NN) attractive and long-range repul¬ 
sive interactions which can stabilize generalized Wigner 



FIG. 1. (color online) (a) Level scheme. Two ground states 
|0) and |1) are coupled coherently through a Raman process. 
The |1) state is weakly dressed to two Rydberg states \Ri) 
and IR 2 ) by a blue- and red-detuned lasers, (b) Interac¬ 
tion potential between two atoms in the Rydberg dressed |1) 
state. Varying the laser parameters and the lattice spacing 
a, we achieve a competing interaction where the NN interac¬ 
tion is attractive and the longer range interactions are repul¬ 
sive. The blue-detuned laser induces the two-atom resonant 
excitation at a distance Rres (dashed line in the figure), (c) 
Tuning the competing interaction by changing the laser de¬ 
tuning A 2 . In the example, two Rydberg states |Ri) = \60S) 
and IR 2 ) = 170*5') of Rubidium atoms are considered. The 
dispersion constants of the respective Rydberg states are 
= 140.4 GHz/xm® and = 882.3 GHz/xm®. Laser 
parameters are given in the following. In (b), we assume 
A 2 = 250 MHz. Other parameters used in plotting (b) and 
(c) are Qi = 30 MHz, Ai = -300 MHz, Q 2 = 22 MHz and 
a = 2 /xm (see text for details). 
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crystals and lead to a new mechanism for the formation 
of a complete devil’s staircase. Here the filling fraction is 
controlled via the NN interaction strength and the emer¬ 
gent devil’s staircase is a union of two sub-staircases — 
a dimer particle staircase and a dimer hole staircase — 
despite of the fact that the microscopic description of 
the system has no explicit particle-hole symmetry. Go¬ 
ing beyond the classical limit we eventually perform a 
qualitative study of the melting of the crystalline states 
through quantum fluctuations introduced by a transverse 
field. In our analysis of the density-density correlations 
we find signatures of a regime featuring quasi long-range 
order before a paramagnetic phase is reached. 

The system .— As shown in Fig. 1(a) we consider atoms 
in two electronic ground states |0) and |1) (typically Zee- 
man or hyperfine states) held in a ID optical lattice. 
They form a pseudo-spin 1/2 particle whose states are 
(Raman) coupled with strength g. One blue-detuned 
and one red-detuned laser are applied simultaneously to 
weakly dress the |1) state with two Rydberg ^'-states \Ri) 
and 1 ^ 2)5 which induces long-range interactions between 
atoms in the Rydberg dressed |1) state [22-25]. The dou¬ 
ble Rydberg dressing in this setting is motivated by re¬ 
cent experimental achievements on Rydberg excitations 
in lattices or microtraps [12, 26-30] and laser dressing of 
ground state atoms to a single Rydberg state [31]. 

In the following, we will elaborate on how to build the 
competing attractive and repulsive interactions in this 
double dressing scheme. The van der Waals (vdW) inter¬ 
action of the Rydberg state \Rj) is /r^ with > 0 
the corresponding dispersion constant {j = 1,2), while 
the inter-state Rydberg interaction can be neglected 
when the two states are far separated energetically [32]. 
The laser (Rabi frequency Clj and detuning Aj) induced 
interaction is given by Uj{r) = /{r^ ± r^) with 

and rj = . Here the 

+ and — sign correspond to the red- and blue-detunted 
laser, respectively. For the dressing by the blue-detuned 
laser, Rres = determines the distance of the two-atom 
resonant excitation when 2Ai + b^i(Rres) = 0 [33, 34] 
[see Fig. 1(b)]. Therefore the resulting interaction I/i(r) 
is attractive for r < Rres and repulsive when r > Rres- 
The dressing of the red-detuned laser generates a barely 
repulsive soft-core interaction U 2 {r). The overall inter¬ 
action potential is given by V{r) = Ui{r) + U 2 {r). As 
depicted in Fig. 1(b), it is divided into an attractive 
(negative) branch and a long-range repulsive vdW part. 
Moreover the relative value W at the attractive branch 
can be tuned by controlling laser parameters, which is 
illustrated in Fig. 1(c). 

In a ID lattice, the number of lattice sites (red cir¬ 
cles) on the attractive branch depends on the length 
Rres- III this work we focus on the case where only the 
NN interaction is attractive. We parameterize the NN 
interaction through W and the long-range tail approxi¬ 


mately through V{k) = CeKak)^ for each k > 2. Here 
Cq « and a is the lattice spacing. The above 

system is described by the Hamiltonian 

H = gJ2^i+Wj2nini+i+ ^ H(|z-j|)(1) 

i i 

where hi = |l)i(l| and s^ = (|l)i(0| + |0)i(l|). Note, that 
we set a as the length unit and V (2) as the energy unit 
without loss of generality. 

Emergent deviVs stairease .—Let us now determine the 
ground state configuration of the model Hamiltonian (1) 
in the classical limit (^ = 0). One notable feature is that 
the crystal formed in the ground state depends on the 
sign of the NN interaction W. In the case of positive W 
the introduction of an excitation is energetically unfavor¬ 
able and the ground state of Eq. (1) is the empty state 
(niG)|0)i). On the other hand in the limit of infinite 
negative W the ground state is the ferromagnetic state 
(niG)|l) 2 ). We are interested in how the ground state 
interpolates between these two cases. 

We first perform a numerical analysis of the infinite 
system. To this end we consider configurations which 
are periodic with period g, i.e. considering rational fill¬ 
ing fractions / = p/g, with p < g, and we set the maximal 
g = 20. For each value of the interaction W we identify 
the ground state numerically as the configuration with 
minimum energy density, i.e., energy of a single period 
divided by its length. Using this method, we obtain a 
staircase pattern connecting the empty state (/ = 0) 
at kF = 0 with the completely filled state (/ = 1) at 
W —1.48, depicted in Fig. 2(a). The data show that 
we never observe isolated excitations, i.e. an atom in 
state |1) surrounded by atoms in state |0), but rather 
polymers of neighboring excitations. This is because ag¬ 
gregating at least two excitations is energetically favor¬ 
able for negative W. We find that the excitations arrange 
in dimers for filling fractions / < 1/2, which is in agree¬ 
ment with the results of Ref. [35]. Once the threshold 
of / = 1/2 is reached the ground state configuration be¬ 
comes • • • 11001100 • • •. In the region of / > 1/2 larger 
polymers form and the numerical results suggest that the 
ground state exhibits a rather intricate crystalline struc¬ 
ture. In the following, we will aim at understanding the 
staircase analytically in both the regimes of / < 1/2 and 
/ > 1 / 2 . 

In the regime of / < 1 /2 — inspired by the numerical 
observations that two excitations are always on neigh¬ 
bouring sites — we formulate an effective dimer particle 
theory to describe the resulting crystalline structure. We 
consider the presence of a dimer particle (DP) on sites 
1) as an effective single particle with associated 
dimer particle number operator p^ = The effec¬ 

tive Hamiltonian of the dimers then reads 

r)pipi+r + WY^pi, (2) 

i r>2 i 
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where W assumes the role of a chemical potential for the 
dimers. Here Vp(r) = 2V{r) V{r 1) V{r — 1) is the 
interaction potential between two dimers separated by r 
sites, which is convex. Note, that for / < 1/2 the convex¬ 
ity of the potential ensures that no neighbouring dimers 
are allowed. We can now employ a method used by Hub¬ 
bard [3] to identify stability regions of different densities 
of dimers with respect to W. The result is plotted in Fig. 
2(b). One observes a remarkable agreement between the 
dimer model and the numerical result by solving the full 
Hamiltonian. 

In the regime of / > 1/2 — unlike the traditional 
staircase where the structure of the staircase in this hole 
sector is trivial, a detailed study must be carried out since 
there is no explicit particle-hole symmetry in our system 
(1). From the numerical data we find that in this regime 
excitations form polymers constituted of R neighboring 
excitations (so-called R-meis). We find regions in which 
only one kind of R-mei is stable, and regions in which 
different kinds of R-meis coexist. The former correspond 



w w 


FIG. 2. (color online) Emergent devil’s staircase, (a) Stair¬ 
case pattern of the filling fraction / versus NN interaction 
strength W. The particle part (/ < 1/2) and the hole part 
of the staircase (/ > 1/2) are indicated by the brown and 
blue shaded areas and merge at filling fraction / = 1/2. (b) 
Magnified view of the dimer particle staircase. The brown 
curve is obtained with the effective dimer particle theory of 
Eq. (2). The black curve (in all panels) is the numerically 
calculated staircase taking into account any rational filling 
fractions f = p/q with ^ < 20 where the range of the inter¬ 
action has been cut off at /c < 1000 lattice sites, (c) At large 
filling fraction (/ > 1/2) the main plateaus (black) corre¬ 
spond to arrangements of a single kind of R-mers. The small 
intermediate plateaus (green) contain mixtures of R-mers and 
R -h 1-mers. The red stars correspond to the predicted tran¬ 
sition points given by Eq. (4). (d) A unified picture for the 
staircase structure at large fillings is obtained by employing 
an effective dimer hole theory of Eq. (5) which agrees excel¬ 
lently with the numerical data. 


to the wider plateaus in the staircase [black data points 
in Fig. 2(c)]. The latter — “mixed” cases — correspond 
to the narrow regions interpolating between two wider 
plateaus [green data points in Fig. 2(c)]. 

As a first approximation, we can construct an R-mei 
theory to describe the staircase as the dominant plateaus 
correspond to different kinds of R-meis (without mix¬ 
tures). Our aim is to determine when the crystal changes 
from R-mei to R 1-mer as we change VF, which is 
achieved by examining the stability of R-meis for a given 
W. Considering an R-mei as a single particle, we define 
the effective filling fraction of R-meis as /eff = P/Q, 
which is related to the real filling fraction / = pjq 
through /eff = p/[R{q — p) + p]- In the following we 
focus on filling fractions with /eff = 1/Q. The energy of 
an P-mer in the thermodynamic limit reads 

oo 

= + (3) 

1 = 1 

where = {R — 1)W + — 0^(0 bond 

energy of the P-mer, which is associated to the mu¬ 
tual interactions of its internal excitations. = 

RV{R—l) [V{x E 1) EV{x — /)] is the interac¬ 
tion energy between two P-mers separated hy x {x = lu), 
where u is the length of the unit cell. 

When the crystal changes from P-mer to P+1-mer, the 
corresponding unit cell length changes from u = QeR—1 
(for P-mer) to u' = Q' R (for P -1- 1-mer). The infinite 
chain can be divided in equal periods on which either 
u' P-mers or u RE 1-mers are disposed, so that we can 
compute the transition point by solving 

uE^^l, = u'E^^K (4) 

We observe numerically that in the regions in which only 
one kind of P-mer is stable the P-mers are homoge¬ 
neously distributed and separated by two holes, giving 
/eff = 1/3 and therefore Q = Q' = 3. This numerical 
observation is confirmed analytically in [36], where we 
show that any other filling fraction makes the resulting 
P + 1-mer arrangement unstable. The comparison of W 
obtained from the simple theoretical model Eq. (4) with 
the numerical data is shown in Fig. 2(c). One can see 
that the transition values of W are in excellent agree¬ 
ment. 

The presence of two neighboring holes in each unit cell 
stimulates us to search for an effective theory of dimers 
of holes (DH). In analogy to the dimers of particles Eq. 
(2), the effective Hamiltonian reads 

jj-DH ^ y^y^\;^(^r)hihi+r - (5) 

i r>2 i 

where Vh(r) = Vp(r) is the interaction energy between 
dimer holes, hi = {1 — hi){1 — n^+i) is the dimer hole 
number operator and ph = 3W E 
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regarded as a chemical potential. We can again study 
the stability regions with a method used by Hubbard 
[3] as the interaction is convex. In Fig. 2(d) we report 
the filling fraction of excitations found in this way, and 
we compare them with the numerical results of the full 
Hamiltonian, finding a perfect agreement in the region 
/ > 1/2. Note by setting the chemical potential = 0, 
no holes are present, and a single polymer extends over 
the whole system corresponding to the aforementioned 
ferromagnetic state. We can then identify the transi¬ 
tion value of W to the ferromagnetic state exactly, i.e., 
Wc = —^0) ~ —1.47994. Note that the emer¬ 
gent staircase is the union of two sub-staircases and is 
complete due to the fact that the long-range interactions 
Vh(r) and Vp(r) in the respective effective models are 
convex [3]. 



FIG. 3. (color online) Quantum melting of the dimer crystal 
with increasing transverse field strength g. Left panel shows 
the transverse magnetization as function of W and 

g, where the simulation was done for L = 60 lattice sites with 
open boundary conditions and the range of the interactions 
has been truncated to 10 lattice sites. The four largest sta¬ 
bility regions f — l,/ = 3/5, / = ll2, f = 2lb are marked by 
dashed lines, see also Fig. 2. Right panel shows the iDMRG 
calculations of the density-density correlation function g^‘^\k) 
vs. k{> 1) for three different values oi g — (0.05,0.15,0.3) at 
W = —0.91 corresponding to the points A,B,G in the left 
panel (see text for details). 

Melting of the stairease .— In the case of a long-range 
repulsive interaction, it was shown in [8] that the intro¬ 
duction of a quantum term leads to a two-stage melting 
process when the excitations are dilute. For weak trans¬ 
verse fields the crystalline structure melts into a float¬ 
ing solid which will eventually melt into a paramagnet 
for higher values of g. Here we explore the qualitative 
features of the melting process in the context of compet¬ 
ing attractive and repulsive interactions for a relatively 
dense crystal. We perform numerical studies by means 
of Matrix Product States implementations of the Density 
Matrix Renormalization Group (DMRG) [41] and its in¬ 
finite version iDMRG [42]. The results are presented in 
Fig. 3. The left panel shows the transverse magnetiza¬ 
tion '^i{sf)/L as function of W and g, where L is the 
length of the chain. The transverse magnetization in¬ 
creases with increasing g. Typically crystals in the mid¬ 


dle of a plateau are more robust against the quantum 
fluctuation, which leads to a series of lobes similar to the 
one found in [8, 43]. 

Here we study the melting of the / = 1/2 lobe which 
has large stability region and is thus naturally easier to 
access experimentally (see [36] for a detailed discussion 
of the melting process for lobes with different filling frac¬ 
tion). We use the density-density correlation function 
g^‘^\k) = {hihij^k) — (which is independent of the 
initial site i for an infinite system) to characterize the 
melting process along the black line in Fig. 3. Perform¬ 
ing a numerical simulation for vanishingly small values 
of g (close to the classical limit) is challenging due to the 
presence of quasi translational invariance which leads to 
a highly degenerate ground state. We consider g = 0.05 
(see point A in Fig. 3), which allows the degeneracy to 
be slightly lifted while still inside the lobe. We find a 
periodic correlation function g^‘^\k) characteristic for a 
crystalline state. In our case {f = 1/2) the ground state 
is quasi-four-fold degenerate and for small values of g has 
the structure \if)g) oc (|..1100..) -f |..0110..) -|- |..0011..) + 
|..1001..))/2 + 0(^^) that allows us to derive the approx¬ 
imate expression g^‘^\k) = cos(/c7r/2) /4 + 0{g‘^) which 
reproduces remarkably well the results shown for point 
A in Fig. 3. 

When g is sufficiently large, the ground state of the 
system is in a paramagnetic phase and the envelope of 
g^‘^\k) shows an exponential decay (point C in Fig. 3). 
By considering intermediate values of g we study the 
melting of the crystalline phase and find that in this 
regime the correlations are significantly longer ranged. 
In the right panel of Fig. 3, we show a typical result 
with bond dimensions of 100 at g = 0.15 (point B in 
Fig. 3). The correlations in the vicinity of the point B 
show similar behavior suggesting that there is a region 
with quasi long-range order extending from 0.1 to 
^:^ 0 . 2 . 

Summary and Outlook .—In conclusion, we have shown 
that a long-range interaction with competing attractive 
and repulsive parts can result in a new mechanism to 
form a complete devil’s staircase. In particular we have 
considered the case of attractive NN interactions which 
lead to a dimer staircase without manifest particle-hole 
symmetry. This suggests an interesting way to control 
the structure of the hole sector of the staircase by tai¬ 
loring the repulsive tail. For example, one will obtain a 
trimer-hole staircase through replacing the longer ranged 
interaction part from the vdW shape (1/r^) to a dipolar 
interaction potential (1/r^). Moreover, since in the out¬ 
lined double Rydberg dressing scheme the range of the 
attractive branch of the potential can be tuned freely, 
and therefore be extended over more lattice spacings, we 
expect that more exotic polymer staircases can be ex¬ 
plored. 
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SUPPLEMENTARY MATERIAL 
I: LIFETIME OF THE DRESSED GROUND STATE ATOMS 

Here we briefly discuss the lifetime of the Rydberg dressed ground state atoms [22] . Lifetimes in the Rydberg states 
|605') and |705') are 252 fis and 410 fis [37], correspondingly. The lifetime in the dressed state |1) is amplified by a 
factor of (2A 1 ^ 27 ^ 1 , 2 )^ — 400 due to the weak mixing with the Rydberg states. With the parameters in the caption 
of Fig. 1 in the main text, we find the effective lifetime is about 0.068 second, which is far longer than other time 
scales in the Hamiltonian. 


II: ANALYTICAL CALCULATION OF THE TRANSITION POINT BETWEEN TWO R-MERS 


In Eq.(4) in the main text we studied the transition between a gas of R-mers and R+ 1-mers assuming the effective 
filling fraction /eff = 1/3 for both gases. Here we want to justify analytically that assumption. As it stands Eq.(4) 
does not constrain the effective filling fractions of the two gases, which can in principle be different. The crucial 
constraint comes from the requirement that at the transition point both the filling fractions of the gas of R- and 
R + 1-mers become unstable. Eor the gas of R-mers, the boundaries of the stability region of a certain filling fraction 
can be determined from the energetics: we call the points at which a gas with one more/less R-mer becomes 
energetically favorable. An analytic expression for /eff = 1/Q was found in [3, 38, 39] and it reads 
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1 

R-1 


^b,int + 

1 = 1 
00 

+ i)E^^\n) - nE^^\ri + 1 ) 


1=1 


(6) 

( 7 ) 


where — l)Vi is the repulsive interaction energy between the constituents of a single R-mer and 

ri = 1{R Q — 1) is the distance between first spins of l-th nearest R-mers. The expressions in Eqs.(6,7) are derived 
under the assumption that the interaction between the R-mers is repulsive, i.e. in our case it breaks down when the 
configuration of the gas can allow two R-mers sitting next to each other, which happens for high filling fractions. We 
need then to restrict to the case ri > R + 1 as for ri = R + 1 the distances ri — 1 = R in Eq.( 6 ) would correspond to 
adjacent R-mers. We plot the devil’s staircase for different R-mers in Eig. 4. 

With the stability region at hand one can find the transition points between different kinds of R-mers as follows: 
firstly we find the regions in W over which two different R-mers are stable, and extract the effective filling fractions 
of the two gases. Second we plug the extracted filling fractions in Eq.(4) to find the transition point. This point 
corresponds to a real transition only if it falls in the overlap of the region of stability of the two gases. 

To illustrate this with an example, let us consider e.g. a filling fraction / = 0.5 which admits both dimers (see 
Eig. 4, black line, /eff = 1/3) and trimers (red line, /eff = 1/4). The condition Eq.(4) then yields W = —1.6858, 
which is inconsistent with the given stability region W G [—0.8319, —0.5759] of trimers. In this way one can show 
that as W is decreased, the chain is filled with dimers only down to W = —1.3066, where the first transition to 
trimers occurs (orange cross at / = 0.5 in Eig. 4). Eollowing this logic it is then straightforward to show, that all the 
remaining transitions between R-mers and RE 1-mers occur at /eff = 1/3 and quickly approach the asymptotic value 
W{R ^ 00 ) = —1.4799. The values of W for some of the transitions are shown in Table I and plotted in Eig. 2c in 
the main manuscript. 


HI: QUANTUM MELTING OF THE / = 2/5 LOBE 

In the main text, we have focused on the quantum melting of the / = 1/2 lobe, which is the largest lobe between 
/ = 0 and / = 1. However, the repulsive van der Waals tail produces an intricate staircase-like structure of the 
crystalline ground states (see Fig.l of the main text). To understand the effects of the repulsive tail, as a first 
approximation, we could keep only the nearest-neighbour attraction and the next-nearest-neighbour repulsion, and 
as such the model reduces to the well-known ANNNI (axial next-nearest-neighbour Ising) model [40], for which only 
three possible ground states exist at ^ = 0, i.e., the / = 0, / = l/2 and / = 1 phases. The transition from the phase of 
/ = 0 to that of / = 1/2 happens at VF = 0 as explained in the main text, and the transition from the / = 1/2 phase 
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FIG. 4. (color online) The devil’s staircase for a gas of i^-mers. The black, red, green and blue staircases correspond to the gas 
of dimers, trimers, pentamers and octamers respectively. The orange crosses indicate the transition points between different 
polymers (e.g. the transition point between dimers and trimers is represented by the orange cross at hlling fraction 0.5). As R 
is increasing the transition value W quickly approaches the asymptotic value W{R = oo) —1.48, which is indicated by the 
orange dashed line. The parts of the staircase for higher filling fractions (after the large plateaus as W is decreased) are not 
shown as they are not captured by our simple theory (see text for details). 


R/iR+l) 

w 

2/3 

-1.3066 

3/4 

-1.4414 

4/5 

-1.4680 

5/6 

-1.4754 

6/7 

-1.4799 

7/8 

-1.4789 

8/9 

-1.4794 

cxo 

-1.4799 


TABLE I. Values of W where a transition between R/{R-\- l)-mers occurs, determined from Eq.(4) with /eff = 1/3. 


to that of / = 1 happens at W/V{2) = —4/3 [40]. Now if one introduces the seemingly small third nearest-neighbour 
repulsion h^(3), which is equal to V(2)/64 due to the 1/r^ law, two additional new phases appear at filling fractions 
/ = 2/5 and / = 3/5. The filling fraction / = 2/5 is stabilized for W G [—0.518, —0.116] which is roughly half of the 
width of the stability region of / = 1/2 {W G [—1.306, —0.522]). Importantly, the filling fraction / = 2/5 is stabilized 
due to the interaction terms going beyond next-to-nearest neighbour. For interaction with the van der Waals tails, the 
physics in the sector with filling fraction / < 1/2 is dominated by pure dimers and the effective interaction between 
these dimers is given by the effective potential Vp rather than by the van der Walls potential itself. 

In Fig. 5 we show the melting process of the / = 2/5 lobe. The qualitative features are analogous to those for 
/ = 1/2 studied in the main text. Inside the / = 2/5 lobe, the crystalline structure is well preserved as indicated 
by the periodic behaviour of the correlation function (point A in Fig. 5). The ground state wave function at 

small g is given by oc (|..11000..) + |..01100..) + |..00110..) + 00011..) + |..10001..))/\/5 + 0{g‘^)^ from which one 
can calculate g^‘^\k = 1, 2, 3,4, 5 • • •) = (1/25,—4/25,—4/25,1/25, 6/25 •••) which matches perfectly the numerical 
results. Experimentally, one can distinguish the different crystalline orders at small g by measuring the different 
periodic behaviours of g^‘^\k). At intermediate g, one observes the quasi long-range order indicated by the algebraic 
decay of the correlations (point B in Fig. 5). Eventually the crystal melts into a paramagnetic phase with exponentially 
decaying correlations (point C in Fig. 5). 
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FIG. 5. (color online) Quantum melting of the / = 2/5 lobe. The left panel shows the transverse magnetization {sf) jL 
in the vicinity of / = 2/5 (see the caption of Fig. 3 in the main text for further details). The right panel shows the iDMRG 
calculations of the density-density correlation function g^^\k) vs. fc(> 1) for three different values of g = (0.02,0.08,0.15) at 
W = —0.32 corresponding to the points A,B,C in the left panel. 
























